Sculpting the Vortex State of a Spinor BEC 
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We use Raman-detuned laser pulses to achieve spatially-varying control of the amplitude and 
phase of the spinor order parameter of a Bose-Einstein condensate. We present experimental results 
confirming precise radial and azimuthal control of amplitude and phase during the creation of 
vortex- ant ivortex superposition states. 
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Many interesting properties of quantum spin fluids 
arise from the existence of topological structures in the 
order parameter. Vortices are among the most funda- 
mental of these topological structures, and play an im- 
portant role in the physics of superfluid helium, super- 
conductors and neutron stars. Bose-Einstein condensates 
(BEC) of alkali metal atoms provide a unique setting 
in which to study superfluid vortices, and have enabled 
great progress in understanding these states [l| . 

Alkali metal BECs have multiple ground state sub- 
levels, and are described by a multi-component order pa- 
rameter. They can carry angular momentum in single- 
component vortices, or in multi-component structures 
such as coreless vortices [H, Q and skyrmions 0, H[ . Cre- 
ating these structures requires precise control of the am- 
plitudes and phases of the components of the order pa- 
rameter. In this Letter, we demonstrate unprecedented 
spatial control of a spinor BEC, using multi-pulse, multi- 
mode stimulated Raman interactions. These results show 
that it is now possible to deterministically create and 
study families of spin textures which have been predicted 
theoretically, but have yet to be directly observed [g, lZ[ . 





FIG. 1: Experimental configuration. A 87 Rb BEC in the 
|F = 1, tuf = —1) state is acted on by simultaneous, coprop- 
agating pulses of a + (light gray) and o~~ (dark gray) polarized 
laser light, with intensities i+, i_. The fields are near reso- 
nant with the 5 2 Pi/2 excited states, and couple the states 
indicated in the energy level diagram (right). 



Our basic experimental configuration has been de- 
cribed previously [9] . For simplicity, the results presented 
in this Letter focus on a specific interaction configuration 
using the F = 1 ground state manifold of 87 Rb, how- 



ever, we have obtained similar results with the F = 2 
manifold. Essential details of the interaction configu- 
ration are depicted in Fig. [H An untrapped BEC is 
acted on by <r + and <j~ polarized laser beams which 
propagate collinearly, parallel to a small magnetic bias 
field (« 1 Gauss). The beams are detuned near reso- 
nance with the 5 2 Pi/ 2 excited states (D\ transitions), 
and satisfy the two-photon resonance condition for the 
Raman transition between the |1) = \F = 1, mp = — 1) 
and |2) = \F = l,rriF = 1) ground state sublevels. The 
individual excited state sublevels which interact with the 
fields are indicated in the diagram by the labels (a-d). 
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FIG. 2: Vortex coupling and interference using beams of 
different orbital angular momentum. Left: the laser pulse 
sequence and beam modes used to interact with the BEC, 
and the state of the BEC as the sequence progresses (a-c). 
An initial laser pulse pair using Gaussian (G) and Laguerre- 
Gaussian (LGj) beam modes causes the initally spin-polarized 
BEC (a) to evolve into a coreless vortex state (b) . Application 
of another pulse pair using G, LGq 1 beam modes, transforms 
this into a coreless vortex-antivortex superposition state (c). 
Right: Experimental images of the BEC (A-C) corresponding 
to schematics (a-c) . The spin components are separated by a 
magnetic field gradient prior to absorption imaging. 

Complex spin textures can be created by applying spa- 
tially inhomogeneous optical fields to a spinor BEC in 
this configuration. Fig. [2] shows a procedure for creat- 
ing first a coreless vortex, and then a vortex superposi- 
tion state by using Laguerre- Gaussian (LG) optical vor- 
tex beams. The coreless vortex is created by applying 
laser fields to the BEC, with the a~ polarized beam in 
an LGq spatial mode, and the <r + polarized beam in a 
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Gaussian (G) spatial mode. The difference in the op- 
tical field modes results in a transfer of orbital angular 
momentum (OAM) to a ring-shaped region of the BEC, 
which is coupled from the |1, — 1) state to the |1, 1) state. 
The center of this rotating ring in the |1, 1) state is filled 
by atoms remaining in |1, — 1); together these constitute 
a coreless vortex state. If a second laser pulse pair is 
then immediately applied which is identical to the first, 
except for a reversal in handedness of the a~ beam mode 
(to LGq 1 ), coupling occurs between this coreless vortex 
state and one of opposite sign. 

The final result of this two-stage procedure is a co- 
herent superposition of a vortex and anti- vortex [1, Q 
in the |1,1) component of the BEC. When the density 
is measured these interfere to give the pattern seen in 
the schematic (c) and the experimental image (C). The 
experimental results (A-C) are obtained by absorption 
imaging along the axis of the BEC after the components 
have been separated by a transverse magnetic field gra- 
dient. It should be emphasized that without this mag- 
netic, Stern-Gerlach separation, the different components 
remain spatially overlapped as indicated by (a-c). The 
interference pattern visible in (C) can be used to confirm 
the winding number of the vortex states, and is also an 
extremely sensitive indicator of spatial inhomogeneities 
in the evolution of the system. This feature is of central 
importance to the results that follow. 
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FIG. 3: (a) Radial intensity profiles for the G, LG beam 
modes used in the experiment, with the BEC density shown 
for comparison. The curves are Gaussian and LG mode func- 
tions fit to data points from azimuthally averaged lineouts of 
images, (b) Experimental images of vortex-antivortex super- 
positions created in |1, 1) by the procedure in Fig. [2] using 
the beam modes of (a). The spiral distortion is caused by in- 
homogeneous light shifts, and increases with beam intensity. 

One of the primary motivations of this work is to pro- 
duce multicomponent coreless vortices in a BEC, with 



sufficient precision to study their dynamics and decay. 
Some of the rich physics of multi-component vortices has 
been glimpsed in previous BEC experiments jloj |. but 
the varieties of states that could be created were lim- 
ited. Raman-coupling of a spinor BEC allows determin- 
istic creation of a wide variety of topological states, par- 
ticularly within the same hyperfine manifold. Many of 
these families of states may be difficult or impossible to 
generate with previous experimental techniques. For ex- 
ample, the structure in Fig[2jB) corresponds to that of 
a half-quantum, or 'Alice' vortex [7]. 

The subtlest challenge in using stimulated Raman in- 
teractions to control a spinor BEC is dealing with spa- 
tially nonuniform light shifts. These can cause spa- 
tially inhomogeneous evolution of the BEC, making state 
preparation and dynamical studies very difficult. The 
spiral distortion of the interference pattern in Fig. EfC) 
is an example of the consequences of these light shifts, 
which arise because the intensity of the beams is not 
uniform over the interaction region. We show here how 
these effects arise, and how to control and compensate 
for them. 

Fig. [31(a) shows the radial dependence of the intensities 
of the beam modes used in the experiment. The Gaus- 
sian beam intensity can be made uniform over the cross 
section of the BEC by increasing its spot size, however, 
the intensity of the LG beam mode cannot be made uni- 
form, because the field vanishes at the central phase sin- 
gularity. The position-dependent intensity of the beams 
causes both the two-photon Rabi frequency, and the dy- 
namic Stark shifts affecting the states, to be position- 
dependent. These light shifts are also state-dependent, 
and cause the phase of the BEC components to evolve at 
different rates at different locations, as if in the presence 
of an inhomogeneous effective magnetic field created by 
the light fields. 

The spiral distortion in Fig. [2](C) is a symptom of ra- 
dial gradients in the relative phase of the BEC compo- 
nents, caused by these so-called Zeeman light shifts. In- 
creasing the effective temporal pulse area of the Raman 
beams magnifies this effect, resulting in tighter spirals in 
the patterns shown in Fig. [3)3. Reversing the handedness 
of the applied LG modes also reverses the spirals seen in 
the interference pattern. Creating a well-formed coreless 
vortex under these circumstances requires the elimina- 
tion of this radial phase gradient. 

Quantitative prediction of the response of the system 
to the laser fields is essential to achieving the necessary 
experimental control. The light shifts that arise during 
the interaction of laser fields with a multi-level atom can 
be complex, however, it is possible to construct a simple 
and accurate model of systems like the one used in these 
experiments. Subject to certain restrictions on detun- 
ing and intensity, one can write a Hamiltonian for the 
system in the interaction picture apply the rotat- 
ing wave approximation, and adiabatically eliminate the 
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excited states [12j. Applying this procedure to the sys- 
tem described by Fig. npresults in the following effective 
two- level Hamiltonian [13], which can be used to describe 
the evolution of the components of the BEC state vec- 
tor in an infinitesimal columnar section of the interaction 
region, parallel to the axis of symmetry. 



H = h 



Xi 1 + 



-xli 



7 ?1 bc 2V / 0+e- i «± 



xih - x\ c I- 



(1) 



In this expression, J + (/_) is the intensity of the <j + 
(<J~) polarized laser field. The two-photon detuning is 
5 = (£"1 — E2)/h— — CO-). E n is the energy of state 
n; cj + (^-) represents the frequency of the <r + (<J~) po- 
larized laser field. The quantity £± is the relative phase 
between the laser fields in the rotating basis defined by 
the frequency difference (cj + — uj-); it is unimportant for 
infinite plane waves, but must be included to describe the 
spatial structure of beam modes carrying different OAM. 
The intensity and relative phase represent the state of the 
optical fields at a point in the plane orthogonal to the 
propagation axis. The atomic external degrees of free- 
dom are represented in this approximation by the phase 
of the spin components; the atomic density is assumed 
to be low enough to neglect mean-field effects. 

The diagonal terms in this Hamiltonian containing the 
parameters \ are the state-dependent light shifts with 
which we are concerned. The off-diagonal term contain- 
ing 77 is the effective two-photon Rabi frequency for the 
coupling between |1) and |2). The sub(super)scripts de- 
note connections to specific ground (excited) states indi- 
cated by labels in Figure [TJ The coefficients determining 
the magnitude and sign of the light shifts (x) and the two- 
photon Rabi frequency (77) depend on the various dipole 
matrix elements, and the detuning of the beams from the 
Zeeman-shifted bare atomic transition frequencies. 

The predictions of this model have allowed us to as- 
sert position-dependent control of the relative phase of 
the components of the BEC, and eliminate the distor- 
tions discussed above. This control is possible because 
the magnitude and sign of the Zeeman light shifts depend 
on the laser detuning. Fig. HJa) shows the calculated rel- 
ative phase shift caused by applying a short pulse of a - 
polarized light to the BEC at different detunings, with 
the intensity and pulse duration held constant. When the 
pulse is detuned between the F' = 1 and F' = 2 states, 
the energy of the |1,— 1) state is decreased relative to 
that of the |1,1). Detuning the field above the F' = 2 
excited state manifold causes the Zeeman light shift to 
reverse sign. In the modified pulse sequence shown in 
Fig. [4f b), an additional pulse of the <j~ polarized field 
is used to create an additional phase shift with identical 
position dependence to that nominally caused during the 
transfer pulses, but with a magnitude and sign that can 
be varied. By choosing an appropriate detuning and in- 
tensity for this control pulse, the inhomogeneity in the 
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FIG. 4: (a) Calculated relative phase shift induced by a a" 
polarized laser field as a function of detuning, for a pulse 
duration of 5 /xs and fixed intensity of 44 mW/cm 2 . (b) Mod- 
ified pulse sequence: A pulse of the a~ polarized field (LG 
mode) is inserted between the 'transfer' pulses, causing an ad- 
ditional spatially inhomogeneous phase shift, with a detuning- 
dependent magnitude and sign. The detunings of the beams 
used in the experiment are indicated by the labels A' and A" 



system can be canceled out, leaving an undistorted core- 
less vortex state. We can confirm this by observing the 
spiral-free interference pattern that would occur in the 
absence of the light shifts. 
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FIG. 5: Experimental images of vortex- ant ivortex superpo- 
sitions in the |1, 1) component of the BEC, created with the 
pulse sequence of Fig.[4jb). (i) A control pulse with the same 
detuning as the transfer pulses tightens the spiral compared 
to the nominal result when no control pulse is applied (ii) . A 
control pulse detuned above the F' — 2 manifold causes the 
light shifts in the sequence to cancel out, resulting in an undis- 
torted interference pattern in the BEC (Hi). Increasing the 
intensity of the control pulse in the 'cancelling' configuration 
overcompensates the radial phase gradients of the transfer 
pulses, creating spirals of opposite handedness (iv,v). 

Fig. [5] shows the experimental results of implementing 
the pulse protocol of Fig. dfb) to increase (i), reverse 
(iv,v), and remove (Hi) the radial gradient in relative 
phase. These are shown in comparison to the nominal re- 
sult generated by the vortex interference pulse sequence 
without the extra control pulse (ii). During the transfer 
pulses, both optical fields are detuned to A'= 0.17 GHz, 
as indicated in Fig. H(a). For the result in Fig. Efi), the 
control pulse detuning is the same as the transfer pulse 
pairs, resulting in an increased spiral in the same direc- 
tion. For results (iii)-(v), the detuning is set to A"= 
0.65 GHz, which results in a Zeeman light shift of op- 
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posite sign, with different magnitudes according to the 
pulse power used. The key result in this set of images is 
(m), in which the spatially inhomogeneous phase gradi- 
ents created in the BEC during the transfer pulses have 
been canceled by the action of the control pulse. 
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FIG. 6: Comparison of experimental vortex interference pat- 
terns and corresponding numerical predictions from a model 
based on the Hamiltonian of Eq. [1] 

The images in Fig. [6] show the experimental results 
for the 'canceling' configuration of Fig. EJ^iii), compared 
against the predictions of a numerical model based on 
Eq.[TJ There is close agreement between the experimental 
and numerical density distributions, which validates the 
approximations used in constructing the model. The nu- 
merical result also provides information about the phase 
of the BEC that is difficult to access experimentally. Note 
that the calculated phase of the |1, 1) state shows the ex- 
pected 7r discontinuity along the density node. The slight 
depression in the phase of both states arises from the un- 
compensated scalar part of the light shifts, which does 
not affect the spin dynamics. This model can be used as 
a guide for creating well-defined perturbations of a core- 
less vortex state, which may prove to be a useful feature 
for future dynamical studies. 

The results summarized in Figs. [5] and [6] demonstrate 
that it is possible to create complex structures such as the 
coreless vortex of Fig. Efb) in a well-controlled manner 
using Raman-coupling methods. The dependence of the 
interaction on excited state structure and field detunings 
allows tremendous flexibility in tailoring the evolution of 
the system for specific experimental goals. The success of 
the approach presented indicates that proposed schemes 
for the storage, manipulation, and retrieval of informa- 



tion in a spinor BEC 14] may be feasible, even in the 
presence of complicated multilevel effects. Furthermore, 
understanding of the principles illustrated in this Let- 
ter are crucial to creating and observing the evolution of 
coreless vortices, for example in an optical dipole trap. 
Further experiments demonstrating the creation of im- 
portant classes of two-dimensional skyrmion states such 
as the Anderson- Toulouse [15] and Mermin-Ho [16] core- 
less vortices are in progress in both the F = 1 and F = 2 
Zeeman manifolds. 
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